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Abstract 

In  this  paper,  a  method  to  monitor  PEM  fuel  cells  internal  temperature  from  surface  measurements  is  presented.  The  aim  of  this  work  is  to 
monitor  fuel  cells  to  prevent  damages  due  to  internal  overheating.  The  measurements  are  taken  at  the  side  of  the  bipolar  plate,  and  heat  flux 
and  temperature  at  the  border  of  the  active  zone  are  estimated.  The  method  is  based  on  sensitivity  analysis  and  inverse  problem  algorithms.  The 
mathematical  formulation  and  algorithm  are  described.  The  model  is  a  transient  heat  conduction  model  in  two  dimensions,  the  inverse  problem  is 
solved  with  an  optimization  method  using  adjoint  equation.  Numerical  test  cases  are  presented  for  graphite  and  steel  bipolar  plates.  The  results 
show  that  internal  temperature  can  be  correctly  estimated.  The  response  time  of  the  method  is  limited  by  the  heat  transfer  rate  in  the  material. 
Therefore,  the  method  is  particularly  appropriate  to  fuel  cells  made  of  steel  bipolar  plates. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Internal  temperature  is  one  of  the  key  parameters  to  monitor 
in  proton  exchange  membrane  fuel  cells.  The  first  reason  is  that 
the  operating  point  efficiency  depends  on  its  value,  the  second 
is  that  a  sharp  raise  in  temperature  can  be  an  indication  that 
direct  combustion  between  air  and  hydrogen  happens  in  the  fuel 
cell.  This  combustion  may  cause  important  damages.  Therefore, 
in  this  paper,  a  method  based  on  inverse  problem  methodology 
to  estimate  internal  temperatures  from  surface  measurements  is 
presented.  The  aim  of  the  study  is  to  investigate  the  capability 
of  monitoring  fuel  cells  with  non-invasive  temperature  sensors. 
Inverse  heat  conduction  problem  is  a  topic  that  has  been  discus¬ 
sed  by  several  authors  [1-4].  In  the  field  of  fuel  cells,  only  a 
few  articles  are  published.  In  references  [5,6]  studies  on  steady 
state  temperature  distribution  at  the  interface  between  the  carbon 
plate  and  the  membrane  electrode  assembly  from  measurements 
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on  the  outer  surface  of  the  end  plate  are  reported.  The  work 
presented  here  deals  with  the  estimation  of  fuel  cell  transient 
internal  temperature  from  surface  temperature  measurements 
taken  at  the  side  of  bipolar  plates.  First,  a  mathematical  method 
based  on  sensitivity  analysis  to  estimate  the  possibilities  of  the 
inverse  problem  methodology  is  described,  then  the  method  to 
solve  the  inverse  problem  of  2D  transient  heat  conduction  is 
presented.  The  methodology  relies  on  approaches  developed  in 
[1-3].  The  paper  is  divided  into  three  sections.  In  the  first  sec¬ 
tion,  the  heat  conduction  model  in  the  fuel  cell  is  formulated. 
The  second  section  describes  the  mathematical  formulation  and 
method,  the  third  part  presents  numerical  results  and  summarizes 
the  conclusions. 

2.  Model  and  formulation 

2.7.  General  description 

A  proton  exchange  membrane  (PEM)  fuel  cell  is  considered. 
Surface  temperature  measurements  are  assumed  to  be  possible 
at  the  side  of  the  bipolar  plates.  The  aim  of  this  work  is  to  esti- 
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Nomenclature 

D  descent  direction 

G  gradient 

h  height  of  area  under  study  (m) 

i  iteration  number 

k  thermal  conductivity  (W  m-1  K-1) 

/  length  of  the  area  under  study  (m) 

L  objective  function 

S  objective  function  including  the  constraints 

t  time  variable 

tf  final  time 

T  temperature  (°C) 

To  initial  temperature  (°C) 

u  heat  flux  (W  m-2) 

v  horizontal  coordinate 

y  vertical  coordinate 

Greek  symbols 

a  thermal  diffusivity  (m2  s-1) 

ft  conjugate  direction  coefficient 

y  descent  step 

8  Dirac  function 

AT  temperature  increment 

A  u  heat  flux  increment 

r  reverse  time  variable 

\jf  Lagrange  multiplier 


mate  the  internal  temperatures  of  the  plate  in  the  active  zone 
of  the  fuel  cell  (Fig.  1).  The  physical  problem  is  formulated 
as  a  heat  transfer  problem.  In  the  border  zone  of  the  bipolar 
plate,  an  unknown  heat  flux  is  considered  to  be  flowing  from 
the  inside  of  the  plate  towards  the  outside.  The  purpose  of  this 
study  is  to  estimate  this  heat  flux  and  the  internal  temperature 
in  order  to  monitor  the  fuel  cell  and  prevent  a  destructive  heat 
elevation. 

2.2.  Problem  formulation 
2.2.1.  Physical  problem 

The  problem  is  a  problem  of  transient  heat  conduction  in  two 
dimensions  with  no  heat  generation  (Fig.  2). 


Conductive  heat  transfer 


Unknown  heat  flux 
and  temperature 

y 1  k 

X 
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temperature 

measurements 


Fig.  2.  Physical  description. 


The  area  under  study  is  limited  by: 

0  <x<l 
0  <y<h 

It  is  considered  that  except  on  the  left  limit,  the  zone  is 
thermally  insulated. 


2.2.2.  Direct  problem  formulation 

With  the  preceding  assumptions,  the  direct  problem  formu¬ 
lation  is  [7]: 

32T  32T  _  1  3T 
3x2  3 y2  a  3 1 


3T 

—k —  =  u(x,  y,  t ),  v  =  0 
3x 
3T 

k —  =0,  x  =  l 
3x 

3T 

—k —  =  0,  y  =  0 
ay 

3T 

k —  =0,  y  —  h 

3 y 

T(x,  y,  0)  =  T0 


(1) 


where  a  denotes  diffusivity,  k  conductivity,  u  unknown  heat  flux, 
Tq  initial  temperature  that  is  supposed  constant. 
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Fig.  1.  Problem  geometry. 


2.2.3.  Sensitivity  problem  formulation 

In  order  to  be  able  to  estimate  what  happens  on  the  left  side 
of  the  area,  we  have  to  check  that  the  heat  transfer  produces  a 
significant  temperature  elevation  on  the  right  side.  This  tempe¬ 
rature  evolution  has  to  be  measurable  by  the  sensors  used.  If  this 
is  not  the  case,  the  estimation  cannot  be  performed.  This  leads 
to  sensitivity  analysis  [1]. 

The  unknown  heat  flux  u  is  perturbed  by  a  small  pertur¬ 
bation  A u  in  (1).  Consequently,  the  temperature  is  perturbed 
by  AT. 
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(1)  is  subtracted  to  the  perturbed  problem.  This  yields  to  the 
sensitivity  problem: 


32A  T  d2A  T 

+ 


1  dAT 


dx2 


dy1 


a  3 1 


dAT 

-k -  =  Au(x,  y,  t ),  x  =  0 

dx 
dAT 


0, 


/ 


dx 

dAT  ^ 

-k -  =0,  y  =  0 

dy 

dAT 

k - =0, 

dy 

A  T(x,  y,0)  =  0 


y  —  h 


(2) 


With  the  results  of  the  sensitivity  analysis,  it  is  possible  to 
design  an  experiment  that  leads  to  a  good  inverse  problem 
solution. 


2.2.4.  Inverse  problem  formulation 

To  estimate  internal  temperatures,  computation  of  the  heat 
flux  at  the  plate  left  side  is  necessary.  The  internal  temperature 
in  x  =  0  is  then  computed  by  solving  the  direct  problem.  Thus, 
the  inverse  problem  of  estimating  internal  temperature  is  divided 
into  two  steps: 

Step  1 :  Compute  u( 0,  y,  t)  from  T(l ,  y,  t) 

Step  2:  Compute  T( 0,  y,  t)  from  u( 0,  y,  t) 

3.  Numerical  computations 

3.1.  Direct  and  sensitivity  problems 

The  direct  and  sensitivity  problem  are  classic  heat  transfer 
problems.  For  numerical  heat  transfer  computations,  we  use  a 
finite  difference  scheme  in  two  dimensions  [8]. 


points,  8  Dirac  delta  function  and  ymeas  coordinates  of  the  mea¬ 
surements.  The  inverse  problem  is  solved  by  minimizing  this 
objective  function  under  the  constraint  that  the  unknown  func¬ 
tion  u(x ,  y,  t)  verifies  the  direct  problem  (1). 


3.2.2.  Minimization  algorithm 

The  conjugate  gradient  method  [9]  whose  algorithm  is  as 
follows  is  employed: 


Do  =  -Go 

Yi- 1  =  Argmin  J(Ui-\  -  yA*-i) 

Ui  =  Ui- 1  +  Yi-iDi-i 

_  {Gf ,  Gj  —  Gj-i) 

Pl~l  IIGi-iH2 

D[  =  —Gi  +  fii-iDi-i 

where  U  denotes  the  unknown  function,  G  gradient,  D  conjugate 
direction,  y  step  size,  (,  )  scalar  product,  and  subscript  i  is  the 
iteration  number. 


3.2.3.  Adjoint  problem- gradient 

The  adjoint  problem  is  derived  by  means  of  writing  the  mini¬ 
mization  problem  as  an  unconstrained  one: 


L{u)  —  f  f  f  (TcompUted(^?  x,  y,  t )  Tmeasured(x,  y,  t )) 

z  J o  Jo  Jo 

x  8(x  -  l)8(y  -  ymeas)dMxdy 


+ 


iKx,  y,  t) 


fd2T  d2T 
\  3x2  3y2 


1  dT\ 
a  dt  ) 


dxdyd^ 


where  j/  denotes  a  Lagrange  multiplier  that  is  solution  of  the 
adjoint  problem. 

The  necessary  conditions  of  stationarity  for  this  functional 
lead  to  the  following  adjoint  problem  [3]: 


3.2.  Inverse  problem 

The  inverse  heat  conduction  problem  (step  2)  is  iteratively 
solved  by  a  gradient  method.  The  gradient  is  computed  using 
an  adjoint  problem.  Among  the  various  methods  to  derive  the 
adjoint  problem,  the  Lagrange  multipliers  are  employed.  The 
following  paragraphs  describe  the  method  used.  More  details 
can  be  found  in  [1-3]. 

3.2.1.  Objective  function 

The  objective  function  is  the  quadratic  residual  between  the 
computed  temperatures  and  the  measured  temperatures: 

i  rh  rl  rtf  2 

S(u)  =  ~Z  /  /  (7eomputed(M!  X,  y,  t )  T'measuredC**  y,  t )) 

JO  JO  JO 

X  S(X  -  l)S(y  -  .Vmeas)tl/dA-d\! 

where  rmeasured  denotes  measured  temperatures  at  the  right 
side  of  the  area,  Tcornputeci  computed  temperatures  at  the  same 


l  djf 
a  dt 


d2\l/  d2\l/ 

ffy2  ^comPuted^’  ^  0 


ed(^  0)^(y  ymeas)^(^)  —  ~ 


d\lf 

-k—  =0,  X  =  0 
3x 
3  \lf 

k—  =  0,  x  =  l 
dx 
d\lf 

-k—  =  0,  y  =  0 
dy 

djr 

k —  =0,  y  —  h 
dy 

\/f(x,  y,  £f)  =  0 


(3) 


By  defining  a  new  time  variable  z  =  tf  —  t,  the  problem  can  be 
integrated. 

The  stationarity  conditions  lead  to  the  gradient: 

G(t)  =  -W,  y,  t) 


Table  1 
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4.  Results 

4.7.  General  description  of  the  numerical  tests 

In  this  section,  the  estimation  of  internal  temperature  of  bipo¬ 
lar  plates  is  numerically  studied.  The  purpose  of  the  numerical 
tests  is  to  determine  if  it  is  possible  to  estimate  an  internal 
temperature  raise  due  to  faulty  conditions  from  surface  tem¬ 
perature  measurements.  Experimental  data  are  simulated.  An 
anomalous  heat  flux  is  supposed  to  happen  inside  the  fuel 
cell  at  v  =  0.  Solving  the  direct  problem  (1)  leads  to  tempe¬ 
rature  values  at  x  =  0  and  at  x  =  l.  Temperature  values  at  x  =  l 
are  considered  to  be  surface  temperatures  issued  from  tempe¬ 
rature  sensors.  Temperature  and  heat  flux  data  at  v  =  0  are  to 
be  reconstructed,  thus  they  are  considered  and  denoted  as  the 
“exact  values”. 

Using  these  simulated  surface  measurements  at  x-l , 
internal  heat  flux  and  temperature  at  x  =  0  are  computed 
using  the  methodology  described  in  Section  3.2.  The  results 
are  referred  as  “computed  values”  in  the  following  para¬ 
graphs. 

The  steel  and  graphite  plates  characteristics  are  displayed  in 
Table  1 .  The  anomalous  heat  flux  applied  at  v  =  0  is  le4  W  m-2. 
Nominal  temperature  operation  is  80  °C.  Thermocouples  K  of 
0.3  °C  accuracy  are  assumed  to  be  used.  Therefore,  the  measu¬ 
rement  of  a  2  °C-elevation  is  possible.  Spatial  discretization  of 
the  numerical  scheme  is  5  x  5. 

4.2.  Numerical  test  case  1:  graphite  bipolar  plates 

In  this  paragraph,  estimation  of  internal  temperature  of  a 
graphite  plate  is  presented. 

4.2.7.  Sensitivity  analysis 

The  influence  of  a  flux  variation  at  the  point  (O,  hi 2) 
on  the  temperature  at  the  point  (/,  h/2)  is  computed  using 
(2). 

A  numerical  test  case  with  50  steps  of  5  s  each  is  conside¬ 
red.  The  results  are  plotted  in  Fig.  3.  They  show  that  a  flux 
variation  at  the  first  time  step  produces  a  sensitivity  coefficient 
that  is  maximal  at  the  20th  time  step.  It  means  that  in  order 
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Fig.  3.  Sensitivity  coefficients  for  50  time  steps  of  5  s,  graphite  plate. 
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Fig.  4.  Surface  temperature  (simulated  measurements)  for  graphite  bipolar  Fig.  6.  Exact  and  estimated  internal  temperature  for  graphite  bipolar  plates, 
plates. 


to  estimate  the  heat  flux  on  one  side  of  the  plate  at  time  step 
1,  the  best  measurement  on  the  other  side  is  the  temperature 
at  time  step  20.  For  the  next  heat  flux  value  (time  step  20), 
the  best  measurement  is  the  temperature  at  time  step  40,  and 
so  forth.  The  sensitivity  coefficient  of  the  40th  time  step  is 
smaller.  The  heat  flux  value  at  this  time  step  cannot  be  well 
estimated  because  the  heat  does  not  have  the  time  to  reach  the 
temperature  sensor.  Of  course,  the  heat  flux  at  the  last  time 
step  cannot  be  estimated  at  all,  the  sensitivity  coefficient  is 
zero. 

So,  this  analysis  shows  that  maximal  sensitivity  is  achieved 
with  a  response  time  of  about  100s.  This  duration  is  too  long 
to  perform  an  efficient  on-line  monitoring  of  the  fuel  cell.  If 
shorter  duration  is  used,  the  temperature  estimation  is  bounded 
to  be  less  precise. 

4.2.2.  Inverse  problem  solution 

In  this  paragraph,  the  inverse  problem  of  estimating  unknown 
heat  flux  and  internal  temperature  values  is  solved  by  means  of 
the  methodology  described  in  Section  3.2.  The  exact  anomalous 
heat  flux  produces  a  temperature  elevation  of  about  50  °C  inside 
the  plate  and  3  °C  at  the  surface  in  50  s  (Fig.  4).  This  simulated 
temperature  surface  is  used  to  compute  the  internal  heat  flux  and 


temperature.  The  computation  lasts  less  than  2  s  on  an  ordinary 
laptop  computer. 

Fig.  5  presents  the  exact  and  computed  heat  flux  values  at 
the  point  x  =  0,  y-h! 2.  The  computed  values  are  found  to  be 
accurately  reconstructed  in  the  first  20  s  of  the  measurement 
interval.  The  relative  error  is  less  than  15%.  After  that,  while 
the  exact  heat  flux  remains  constant,  the  computed  heat  flux 
strongly  decreases.  Heat  transfer  in  graphite  is  too  slow  to  pro¬ 
duce  significant  surface  temperature  elevation  in  less  than  30  s. 
Therefore,  the  algorithm  lacks  information  and  is  not  able  to 
reconstruct  heat  flux  values  in  the  last  part  of  the  interval.  The 
accuracy  decreases.  This  result  is  in  accordance  with  the  sensi¬ 
tivity  analysis.  Fig.  6  presents  the  exact  and  computed  internal 
temperature.  As  before,  the  values  are  accurately  computed  in 
the  first  half  of  the  time  interval.  The  maximum  relative  error 
between  the  exact  and  computed  values  remains  under  5%  for 
28  s.  This  accuracy  is  better  than  the  accuracy  obtained  for  the 
heat  flux.  It  is  due  to  the  smoothing  effect  of  the  heat  equation 

[1-4]. 

These  observations  led  to  the  fact  that  it  is  necessary  to  take 
into  account  the  duration  of  the  heat  transfer  process  to  achieve 
accurate  estimation  of  internal  values.  In  this  example,  if  the 
temperature  measurements  interval  is  of  50  s,  the  estimation  of 
the  internal  values  is  accurate  only  for  the  20  first  seconds. 
So,  the  proposed  method  produces  5%  accurate  temperature 


Fig.  7.  Sensitivity  coefficients:  50  time  steps  of  1  s,  steel  plate. 


S.  Begot,  J.M.  Kaujfmann  /  Journal  of  Power  Sources  178  (2008)  316-322 


321 


Time  (s) 

Fig.  8.  Surface  temperature  (simulated  measurements)  for  steel  bipolar  plates. 


Fig.  9.  Exact  and  estimated  heat  flux  for  steel  bipolar  plates. 


values  with  a  50-s  measurement  duration  and  a  30-s  response 
time. 

4.3.  Numerical  test  case  2:  stainless  steel  bipolar  plates 

In  this  paragraph,  the  bipolar  plates  are  now  supposed  to  be 
made  of  stainless  steel.  The  steel  characteristics  are  presented  in 


Fig.  10.  Exact  and  estimated  internal  temperature  for  stainless  steel  bipolar 
plates. 


Table  1 .  As  the  thermal  conductivity  of  steel  is  much  higher,  an 
improvement  in  the  response  time  of  the  temperature  estimation 
is  expected. 

4.3.1.  Sensitivity  analysis 

Fig.  7  presents  the  sensitivity  coefficients  for  steel  plates.  The 
coefficients  are  higher  than  those  for  graphite  plates.  Only  the 
coefficient  of  the  40th  time  step  is  strongly  smaller. 

4.3.2.  Inverse  problem  solution 

The  internal  heat  flux  and  temperature  are  estimated  from  sur¬ 
face  simulated  measurements.  The  exact  heat  flux  is  le4  W  m-2 
as  in  the  preceding  test  case.  As  the  thermophysical  properties 
of  steel  and  graphite  are  different,  temperature  evolution  in  the 
two  cases  differs.  Globally,  the  increase  in  temperature  for  steel 
plates  is  less  important,  but  the  heat  transfer  is  quicker.  Fig.  8  pre¬ 
sents  the  simulated  surface  temperature  of  the  bipolar  plate.  The 
estimated  heat  flux  and  temperature  are  plotted  in  Figs.  9  and  10. 

As  expected,  a  clear  improvement  in  heat  flux  and  tempera¬ 
ture  estimation  is  observed.  For  a  50-s  measurement,  the  heat 
flux  estimation  is  accurate  for  70%  of  the  measurement  interval. 
The  internal  temperature  is  computed  with  an  accuracy  better 
than  5%  for  48  s.  Before  35  s,  the  accuracy  of  the  temperature 
estimation  is  even  better  (2%).  This  result  is  in  accordance  with 
the  sensitivity  analysis.  Heat  transfer  in  the  steel  is  fast  enough 
to  produce  a  significant  surface  temperature  elevation  in  about 
15  s.  Therefore,  the  algorithm  has  sufficient  data  to  produce 
accurate  results.  An  accurate  estimation  of  internal  tempera¬ 
ture  can  be  carried  out  with  a  2%  accuracy  and  a  15-s  response 
time. 

So,  the  proposed  method  enables  us  to  estimate  internal 
temperature  with  a  predicted  accuracy  and  response  time.  The 
performance  of  the  method  is  linked  to  the  thermophysical  pro¬ 
perties  of  the  material,  so  it  is  particularly  appropriate  for  fuel 
cells  made  of  steel  bipolar  plates. 

5.  Conclusion 

In  this  paper,  a  method  to  estimate  internal  fuel  cell  inter¬ 
nal  temperature  from  surface  measurements  is  presented.  The 
aim  is  to  monitor  fuel  cells  to  prevent  damages  due  to  inter¬ 
nal  overheating.  The  measurements  are  taken  on  the  side  of 
the  bipolar  plate,  and  the  heat  flux  and  the  temperature  at 
the  border  of  the  active  zone  are  estimated.  The  method  is 
based  on  sensitivity  analysis  and  inverse  problem  algorithms. 
Numerical  results  show  that  internal  temperature  can  be  cor¬ 
rectly  estimated.  The  response  time  of  the  method  is  limited 
by  the  heat  transfer  rate  in  the  material.  Therefore,  the  method 
is  particularly  appropriate  to  fuel  cells  made  of  steel  bipolar 
plates. 
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